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We theoretically study the propagation of sound waves in GaAs/AlAs superlattices focussing on 
periodic modes in the vicinity of the band gaps. Based on analytical and numerical calculations, 
we show that these modes are the product of a quickly oscillating function times a slowly varying 
envelope function. We carefully study the phase of the envelope function compared to the surface of 
a semi-infinite superlattice. Especially, the dephasing of the superlattice compared to its surface is 
a key parameter. We exhibit two kind of modes: Surface Avoiding and Surface Loving Modes whose 
envelope functions have their minima and respectively maxima in the vicinity of the surface. We 
finally consider the observability of such modes. While Surface avoiding modes have experimentally 
been observed (Phys. Rev. Lett. 97, 1224301 (2006)), we show that Surface Loving Modes are likely 
to be observable and we discuss the achievement of such experiments. The proposed approach could 
be easily transposed to other types of wave propagation in unidimensional semi-infinite periodic 
structures as photonic Bragg mirror. 

PACS numbers: 



I. INTRODUCTION 

Propagation of waves in periodic structures have been 
studied for decades^. Indeed, the periodicity leads to 
band gaps where propagation is forbidden. This prop- 
erty is very general and could be observed in many dif- 
ferent fields. Electrons in crystals experience the peri- 
odic atomic potential leading to electronic band gaps 2 . 
Electromagnetic waves cannot propagate in a spatial pe- 
riodic dielectric constant system if their frequencies fall 
in the band gap: Photonic crystals 3 -^ are based on this 
statement. Elastic waves in solids could also endure a pe- 
riodic Young modulus and/or a periodic density. Studies 
on such systems are currently numerous^i 7 - and lead to 
a new field of physics namely Phononics 8,9 . 

In periodic structures, two kinds of vibrations can be 
pointed out: 1) extended eigenmodes that propagate 
through the system; their frequencies fall out of the gaps 
and they satisfy the Floquet-Bloch theorem; 2) localized 
modes that can be found around isolated defects includ- 
ing surfaces^surface modes belong to this category). 
Due to their localized character, these modes are not so 
sensitive to the periodicity of the underlying Crystal and 
thus their frequency can take any value, in particular 
within the forbidden band gaps. 

Extended acoustical vibrations in superlattices (SL) 
have been described quantitatively in the frame of the 
elastic model in the middle of the century-ii. The su- 
perlattice, with its long period compared to the under- 
lying crystals periods, results in a much shorter recip- 
rocal Brillouin zone and hence a multiple folding of the 
initial acoustical branches. Experimentally, observations 
of these folded acoustical vibrations have first been per- 
formed using Raman scattering experiments for which 
selection rules are now well understoo d 12 ^ 3 . Since the 



nineties, time-domain optical experiments have investi- 
gated folded acoustical vibration o 14 ' 15 ! 16 ' 17 . They exhibit 
two types of modes: phonons with finite wave vectors 
(q =/= 0) and zone-center modes (q = 0) in the reduced 
Brillouin zone. 

Recently, Trigo et al^ studied both theoretically 
and experimentally these zone-center modes in phononic 
semi- infinite GaAs/AlAs SL yielding two conclusions. 
First, modes near the Brillouin zone-center observed by 
pump-probe experiments have not rigorously a zero wave 
vector, but instead an almost zero wave vector q s» 0. 
Second, these modes present a slowly varying envelope 
wave function with an amplitude minimum in the vicinity 
of the surface: they claim these modes "have a tendency 
to avoid the boundaries, irrespective of the boundary 
conditions" and for this reason, refer to them as Surface 
Avoiding Modes (SAM). Note that these surface avoid- 
ing modes have also been encountered near photonic 19 
and acoustic^ band gap materials and their nature is 
very general. In this work, we study in detail acoustical 
modes in the vicinity of Brillouin zone center, focussing 
on their avoiding character near the GaAs/AlAs SL sur- 
face. 

Our work is based on analytical studies corroborated by 
numerical simulations. Analytical calculations are based 
on the analogy between the propagation of waves in a 
periodic structure and the spatial parametric oscillator. 
Since this analogy is not well-known, we devote Sect.HTlto 
prove the fruitfulness of such analogy to describe waves 
in periodic structures in the vicinity of gaps. Considering 
the propagation in an infinite SL of longitudinal modes 
in the vicinity of band gaps, we show that these modes 
are the product of a fast oscillating function and a slowly 
varying envelope function. We derive the dispersion di- 
agram in the vicinity of band gaps using the parametric 
oscillator equations. Then considering the case of a semi- 
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infinite SL in Sect. Hill we give the analytical expression 
of longitudinal waves assuming free boundary conditions 
and considering the dephasing of the SL compared to its 
surface. We show the existence of Surface Avoiding and 
Surface Loving Modes depending if the envelope function 
is minimum or maximum near the SL surface. We then 
discuss in detail the avoiding or loving character of modes 
near the center of the Brillouin zone. In Sect. IIV1 we dis- 
cuss on the electron-phonon coupling of these modes and 
show both Surface Avoiding and Loving Modes are likely 
to be observable using time-domain optical experiments 
or Raman and Brillouin scattering experiments. Finally, 
Sect. W\ will be devoted to a discussion on our results. 



II. PARAMETRIC OSCILLATOR 

A. Parametric oscillator analogy 

We consider the propagation of a sound wave in an 
unidimcnsional infinite SL made of 2 solid materials: re- 
spectively named A and B. cIa and ds are the width of 
the layers, and d — d& + ds is the period of the SL. If u 
is the atomic displacement, in each material sound waves 
are solutions of: 
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where Ci is the sound speed in the materials (i — A, B). 
To determinate sound waves in a SL from Eq. J]}, one 
needs to apply at each interface, the continuities of the 
displacement field and of forces dF — adS acting on each 
elementary interface dS, where a is the stress tensor. 

In a SL, the displacements can be determined exactly 
by the transfer matrix method^. Though very powerful, 
and clearly grounded, this method suffers from its numer- 
ical nature. It is not easy to get analytical expressions 
for finite size superlattices. 

In this work, we will use the analogy between the prop- 
agation of sound waves in periodic media and the para- 
metric oscillator. As it will be shown, this analogy in- 
volves severals approximations we will justify. However, 
it has two benefits. First, it yields analytical expres- 
sions of the displacements in closed form, with relatively 
simple calculations. Second, and especially, it allows to 
point out and to understand the main physical phenom- 
ena near the Bragg reflexions in the SL. We stress that 
the present treatment is complementary to the transfer 
matrix method and does not pretend to challenge it for 
quantitative analysis. 

Finally, one must note that the present method could 
be applied to the propagation of light in photonic crystals 
or electrons in crystal s 22 ! 23 since their wave propagation 
equations are similar. 

From now on, we consider a harmonic, plane, longi- 
tudinal sound wave of angular frequency to: u(z,t) = 



layers. Continuity of forces acting on each interface im- 
poses the continuity of Cn (z) 9w } where Cu is the stiff- 
ness constant, for a longitudinal sound wave propagating 
along z. In the following, we will focus on GaAs/AlAs 
superlattices, because they are the most experimentally 
used superlattices. Stiffness constant contrast between 
GaAs and AlAs is very weak: ~ 0.014 so that the 

continuity of Cu(z) dw ^ at each interface reduces to the 

continuity of a " W . I n such a case, the determination of 
sound waves in the SL can be directly achieved by the 
resolution of Eq. (Q]) considering only the variation of the 
sound speed since u(z) becomes a C 2 function 2 ^. Con- 
sequently, since the sound speed in a GaAs/ AlAs SL is 
a periodic function of z, we may expand as a Fourier 
series £ = a + £"=1 a m cos(^p) + b m sin(^p), so 
that sound waves in GaAs/ AlAs superlattices are solu- 
tions of: 
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the primitive vector of the 
— the reduced 

a 



We now introduce G = =§■ 
reciprocal lattice, c m = Sat a nd d 

amplitudes of the mth harmonic, and fco(w) = ujy/aQ the 
zero order wave number i.e. the wave number a wave 
would have in the absence of any periodic modulation of 
the sound speed. ko(u>) — u^/oq is the usual dispersion 
relation for acoustic phonons in an homogeneous medium 
with mean sound speed cq = ^= . Eq ^ then reads like 
a parametric oscillator equation in spatial coordinates: 
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, cos(mGz) + d m sin(mGz) 



u = 



(3) 



B. Single Harmonic Approximation 

Eq. §3§ includes all harmonics of the inverse square 
sound velocity. In order to understand the effect induced 
by the different harmonics, let us first consider the effect 
of an isolated harmonic. Eq. (J3j> , in this case, reduces to 
a Mathieu equation 24 of the form: 



dz 2 



+ k 2 [1 + a cos(kez)] u = 



(4) 



u(z)e 



along the z-axis, perpendicular to the SL 



k e is the wave number of the periodic excitation. General 
solutions of Eq. ^ are periodic except in some regions 
of the plane [ko/k e , a): they are then the product of an 
oscillating function with a linear combination of an in- 
creasing and a decreasing exponential. These non peri- 
odic solutions correspond to the resonances of the para- 
metric oscillator or band gaps of the SL. Fig. [T] presents 
the stability diagram of Eq. (H|) limited to the first four 
resonances as a function of ko/k e and a: regions of non- 
periodic solutions are dashed. This diagram has been 
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FIG. 1: Phase diagram of the Mathieu equation. Parameters 
inducing non periodic solutions are in dashed regions ("para- 
metric resonances" or "band gaps"). 



calculated numerically using the resolvent of Eq. and 
the Liouville and Floquet-Bloch theorems^. 

Regions of non periodic solutions form in the plane 
(fco/fc e ,a) bands that, when a tends to 0, converge 



to points f 2 - = § with n 6 IN*, i.e. wave numbers 
k e = 2fc , ko, §fco, . . . etc. Moreover, for a <C 1, it can 
be shown^i that the band width of the nth resonance 
is a decreasing function of n proportional to a". As a 
result, the larger resonance is obtained for n = 1 i.e. for 
k e = 2fco, the most famous resonance of the parametric 
oscillator, corresponding to an excitation twice faster 
than the zero order wave number. 



Hence, going back to Eq. ([3]), varying the value of 
ko, each harmonic m of the excitation at wave num- 
ber mG, if treated independently, is expected to create 
band gaps around = ~. Thus, band gaps occur for 

2gi = nm = I G IN* . Fig. [2] schematically presents the 
band gaps independently created by each of the first fifth 
harmonics. Excepted the first band gap (1 = 1) due ex- 
clusively to the first harmonic (m = 1), all other band 
gaps result from at least two harmonics. However, as 
a first approximation, since the larger resonance due to 
the mth harmonic is for n = 1 i.e. around = 5, 
we consider that the Zth band gap of Eq. §5§ is essen- 
tially due to the first resonance (n = 1) of the Ith har- 
monic (m = I). We call this approximation the Single 
Harmonic Approximation (SHA). This approximation is 
equivalent to treating the modulation of the sound speed 
as a first order perturbation: |c m | <C 1 and \d m \ <C 1 
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FIG. 2: Schematic illustration of band gaps produced by each 
harmonic. 



Vm G IN*. Our main motivation for the SHA is the possi- 
bility to derive analytical calculations and the possibility 
to produce good analytical approximations of solutions 
of Eq. ((3]) around a given gap. Of course, the validity 
of SHA will be discussed below in the considered case. 
Note that values of ko for which there are no periodic 
solution, correspond to a range of frequencies u> which is 
the usual band gap of the SL. Indeed, fco has been defined 
by kg — loJoo in Sect. Ill Al 



C. Infinite medium: band gap 



In order to test the validity of SHA, we solve Eq. ([3|) 
near the resonance^ at oj^q defined by fc (w^g) = G 
i.e. the first resonance due to several harmonics (m = 1 



and m — 2). The expression of uj^q reads: u. lsa 

(1) 

The superscript "(2)" in ui B q refers to the second band 
gap. Using the SHA, this band gap is only produced by 
the harmonic at 2G in the Fourier series of \. We thus 
study a reduced equation taking into account only that 
harmonic: 



,(2) 



G 



^— + kl [1 + c 2 cos(2Gz) + d 2 sin(2Gz)l u = (5) 
dz 2 



Where we have simplified the notation of ko(u>) using 
ko- The SHA will be checked a posteriori. Switching 
to complex number representation and introducing p 2 = 
Eq. dU reads: 



c 2 —id2 



-p 2 e i2Gz +p 2 e~ 



i2Gz 



)U = 0, 



(0) 



where p 2 denotes the complex conjugate of p 2 . Similarly 
to the resolution of the parametric oscillator equation by 
Landau and Lifchitz^i, we change the coordinate system 
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to A(z) and B(z) denned by: 
u(z) = A{z)e lGz 



B(z)e 



-iGz 



(7) 



where A(z) and B(z) are slowly varying functions. Intro- 
ducing Eq. in Eq. © and neglecting fast oscillating 
terms (at wavevector 3 G) as well as the second deriva- 
tives of A(z) and B(z), the identification of coefficients 
of plane waves at G and — G reads: 
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whose general solutions are: 

A{z) = A+e iKZ +A-e~ lKZ 
B(z) = B+e lKZ + B_e- tKZ 

where n is given by: 

(k 2 -G 2 ) 2 \ P2 \ 2 kt 



(8) 



(9a) 
(9b) 
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8 characterizes the detuning between the excitation wave 
number G and the zero order wave number ko- 
The system undergoes a bifurcation when k 2 changes its 
sign as a function of to: 

• Either k 2 < 0. Solutions of Eq. ^ are then hyper- 
bolic: A(z) and B(z) are a linear combination of 
increasing and decreasing exponentials. And there- 
fore, no propagation can occur in the SL. The con- 
dition k 2 < defines the band gap of the SL around 

^BG- 

• Or k 2 > 0. Solutions of Eq. ((9|) are sinusoidal and 
the general solution u(z) is thus periodic. Note that 
in this case, the spectrum of u(z) will be composed 
of a mixing of G + k and G — k plane waves. u(z) is 
then the product of a quickly oscillating function at 
G, times a slowly varying functions at k. Moreover, 
Eq. (fT0|) establishes a relation between oj (through 
ko(u>)) and wave vectors G± k i.e. the dispersion 
relation of the SL around the considered gap. 

Note that the band gap only depends on the squared 



magnitude of p2- \p%\ 



ci+di 



i.e. on the amplitude of 



the second harmonics at 2G. 

Besides, we already pointed out the similarities of sound 
waves propagation in SL and electron waves propagating 
in crystals 22 . Usually, when studying the independent 
electrons in a weak crystalline potential, one applies the 
first order perturbation theory to determine the width 
and position of electronic band gap 2 : the expression of 
the band gap width is analogous to Eq. (|9"aj) and could be 
derived using the parametric oscillator analogy and the 
SHA^. 



D. Dispersion diagram: numerical study 

In this part, we will numerically check our analytical 
results of Sect. Ill Cl and the SHA applying our predictions 
on a GaAs/AlAs SL with cLa = d,Q a A s = 5.9 nm and 
<Ib = d-AiAs = 2.35 nm. Such a SL has been used in an 
experimental setup by Trigo at al—L For these materials, 
longitudinal sound speeds are^i: 



ca = c G aAs = 4726 ms 1 
cb = caias = 5630 ms -1 



(12a) 
(12b) 



These parameters will be used in all the numerical exam- 
ples given in this article. 

Fig. [3] compares the dispersion diagram obtained 
from Eq. (fTTj)) (SHA approx.) and the transfer matrix 
method2i(exact solutions) for this SL. Fig.[3]is restricted 
to the two first band gaps respectively due to the har- 
monics at spatial frequencies G and 2G. We would like 
to emphasize that in Eq. (fT0|) . the first band gap is ob- 
tained using coefficient p\ = Cl ~ tdl (instead of P2) and 
the second bang gap using only p2 = C2 ~^ d2 whereas the 
transfer matrix method considers all harmonics. 
According to Fig. [31 our predictions of band gaps with 
the SHA are in good agreement with the exact calcula- 
tion of the transfer matrix method. Moreover, dispersion 
diagrams around the gaps also agree very well. They 
slowly diverge when moving away from the gaps: the A 
and B functions are varying faster and thus their sec- 
ond derivatives are not negligible any more. In addition, 
the development in Eq. ([7]) neglecting harmonics of order 
higher than 2 and keeping only the main term around 
wave number G is no longer appropriate. Finally, the 
SHA, keeping only one harmonic in Eq. ([SJ naturally be- 
comes inadequate when moving too far away from the 
gap. 

Fig. [3] validates our SHA for the first and second band 
gaps; and it is easy to convince oneself that it is also 
true for any of the SL band gaps in a first order ap- 
proximation. Below, we will focus on periodic solutions 
around the second band gap: SL eigenmodes in the vicin- 
ity of this gap will thus be well described by solutions of 
Eq. ©. 



III. SAM AND SLM 



In the previous section, we have studied the general 
solution of Eq. ^ and derived the dispersion diagram 
of the SL. Since experiments always involve a finite SL, 
we now consider the effects of the presence of a free sur- 
face at z = with definite boundary conditions on SL 
eigenmodes. 
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FIG. 3: (color online) Dispersion diagram of a GaAs/AlAs SL 
calculated using the transfer matrix method (dashed black) 
and predictions of the gaps using the parametric oscillator 
analogy and the SHA: (red) first gap calculated from the first 
harmonic of the inverse squared sound speed and (blue) sec- 
ond gap calculated from the second harmonic. We report 
using green horizontal dotted lines the positions of the lih 
band gap uJg G = 2 ^L_ . P arameters corresponding to the SL 
described in section Hi Dl have been used. 



Semi-infinite medium: surface effects 



More precisely, we intend to study periodic solutions of 
Eq. ([3]) in the vicinity of the band gap at oj bg using free 
boundary conditions i.e. u(0) — uq and ^f(O) = 0. In 
the following, we will use the abbreviation NBPM (Near 
Bang gap Propagative Mode) in reference to such modes. 
Figure 2] plots ^ as a function of z. The dephasing r 
of the SL compared to z = is defined on Fig. 0J As 
shown below, t will be a key parameter in our study. 
Coefficients ao, c m and d m in Eq. ([3]) are calculated from 
the Fourier components of \: 
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FIG. 4: (color online) Schema of the SL compared to the 
surface z — 



Where a — r/d is the reduced dephasing and 7 = dA/d 
is the cycle ratio of the SL. The mean sound speed cq is 
defined by Eq. (fHa|) . 

Since we study NBPM in the vicinity of the second 

(2) 

band gap around u BG , it is instructive to consider the 
evolution of P2 as a function of the cycle ratio 7 for dif- 
ferent values of the reduced dephasing a as shown in 
FigEl P2 vanishes twice as the cycle ratio 7 varies from 
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FIG. 5: (color online) Plot in the complex plane of p2 (|13d|l . 
for 7 G [0, 1]. The first case considered is a = (black). The 
curve is not a single circle, since the average sound speed en 
depends on 7, as seen in Eq. (|13a[) . The reduced dephasing a 
of the SL merely amounts to a rotation of the curve around 
0. The dot circles correspond to the SL taken as an example 
in this paper with a cycle ratio 7 = ysL- For a = a% = 
1/2 — jsl/2 or a = Qi + 1/2 = 1 — jsl/2 (red) or a = 02 = 
3/4 - 7sl /2 or a = a 2 + 1/2 = 5/4- 7sl /2 (blue) , this brings 
P2(7sl) on the real axis. 

to 1. So there are two values of 7, independent of a, 

(2) 

for which the SHA predicts no gap opening around oj bg - 
Note that a gap might actually be opened through the 
second parametric resonance n = 2 produced by the first 
harmonic m — 1 (cf. Fig. 

For a given 7, the effect of a is merely to rotate the 
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figure by — 2-Kma : as seen in Eq. (|13d[) . Note that experi- 
mentally, most of the SL have a — and that for a = 0, 
P2 as a function of 7 never cross the real axis except when 
P2 = 0: we will illustrate below to what extend the case 
P2 non null and real is relevant. 

Details of the resolution of Eq. © are given in Ap- 
pendix |XJ Since we are focussing on periodic modes, i.e. 
outside the gap, k 2 in Eq. (TIT))) is positive. Let us intro- 
duce ip the argument of P2 ■ 



P2 = |P2| 



itp 



(14) 



The solution u(z) satisfying the free boundary conditions 
w(0) = mo and ^(0) = reads 

u{z) = 2A a [(1 - C) cos(Gz + VV 2 ) cos(kz + <p) 

-(l + ()sm(Gz + i)/2)sm(Kz + 4>)} (15) 

where £, <f> and Aq are defined by 



c 
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tan<^> = — tan('0/2) 
A = 



G(l-C) + «(l + 



G(l 



■C) + «(i 
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(16) 

(17) 
-418) 



(1-C) cosf cos0- (1 + C)sin- 



u(z) is real and u(z, t) is a stationary wave, as expected 
since the surface is a perfect mirror for acoustic waves. 

Eq. (|15[) is linear combination of two terms each involv- 
ing a quickly oscillating function at G and a slowly one at 
k corresponding respectively to Bloch and envelope wave 
functions. The parameter £ governs the relative contri- 
butions of each term and hence their dominant character 
which we will discuss in Sec. IIIIB 21 The parameter <j) 
is a dephasing of the envelope function with respect to 
the SL surface: it will be discussed in Sect. IIII B 31 Aq 
governs the global amplitude of the displacement field. 
There are in general beatings in the solution u(z), as em- 
phasized by the analytical expression for its root mean 
squared envelope it rms , derived in Appendix IBl 



< ns (z) = 2A£ 1 + C - 2Ccos2(kz + ■ 



(19) 



where the fast oscillating terms at G present in Eq.JT 
have been smoothed out. 



1. Remarkable cases. Pure SAM and pure SLM 

• P2 is real andp2 > (a — 0,2 or 02 + 1/2 on Fig. [5]) . 
Therefore from (Q3]), ip = and from |T7|), <j> = 0[n]. 

— Just above the gap, we have <5 « \p2\ (cf 
Eqs. (TO]) and CLU)), implying ( « l(Eq. (fTH|> ). 
Thus from Eq. (|15p . the main contribution to 
u(z) is proportional to sin(Gz) sin(nz): the 
envelope of the vibration mode is null at the 
surface. Following the denomination of Mer- 
linSii, we call such a mode a Surface Avoiding 
Mode (SAM): these SL eigenmodes shy away 
from the boundaries. 

— Just below the gap, S w — \p2\ and thus 
£ w —1. The main contribution to u(z) is 
cos(Gz) cos(kz): The envelope amplitude is 
maximal at the surface. In view of that, we 
refer to such SL eigenmode as Surface Loving 
Mode (SLM). 

• P2 is real andp2 < (a = ai or ai + 1/2 on Fig. [5]), 
ip = 7r and thus (f> — tt/2[tt]. 

— Just above the gap, S « \P2\- The main contri- 
bution to u(z) is then sin(Gz + 7r/2) sm{nz ± 
7r/2) = ± cos(Gz) cos(kz). Therefore, the SL 
eigenmode corresponds to a Surface Loving 
Mode. 

— just below the gap 5 « — The main con- 
tribution to u(z) is cos(Gz + 7t/2)cos(kz ± 
7r/2) = ± sin(Gz) sin(«:z): Therefore, the SL 
eigenmode corresponds to a Surface Avoiding 
Mode. 

Note that we have chosen to base the latter discussion 
on Eq. (fT5"|) , but Eq. ([IT?]) could also be used and leads to 
the same conclusions. 

If the distinction between SAM and SLM is obvious in 
the cases where P2 is real, in the general case we need a 
more precise definition: we will speak about SAM or SLM 
depending if the envelope amplitude at the surface (z = 
0) is smaller or higher than its value a quarter period later 
at z = 4r = 7T- where X K is the wavelength associated 
with the envelope wave function. With this definition, 
according to Eq. fH]), knowing ( and cf> is enough to 
determine the SAM or SLM character. 



B. Surface Loving and Avoiding Modes 

In this section, we discuss the effect of the reduced de- 
phasing a of the SL compared to the surface z = on 

(2) 

NBPM near the second gap around ^ B q- We especially 
focus on the phase of the envelope wave function at wave 
vector k. Examining the general case requires a numer- 
ical study. However, two limiting cases where P2 is real, 
are remarkable. 



2. Beatings contrast 

Let us focus on £ which, as previously quoted, governs 
the relative amplitude of the two contribution in Eq. (|15p 
and hence drive the beatings contrast, which is defined 
in Eq. (|B4p and plotted in Fig [6] as a function of the fre- 
quency v — From Fig. [S] the contrast reaches its 
maximum value of 1 only at the precise band gap edges. 
In such a case, the envelope function exactly vanishes at 



its minimum: one could then speak about perfect SAM 
or SLM. Moving away from the gap, the contrast tends 
to giving the same weight to both terms in Eq. (fT5|) . 
These beatings are also present in an infinite SL, just like 
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FIG. 6: Beating contrast C of the displacement u, defined in 
Eq. ()B4[l as a function of the frequency v — C equals £ 
above the gap and — below. Parameters corresponding to 
the SL described in section III Dl have been used. Dot circles 
mark a beating contrast of 0.5 for frequencies 0.5807 THz and 
0.6182 THz. 

in a temporal parametric oscillator. Away from the reso- 
nance condition, the excitation goes successively in-phase 
and out-of-phase with the movement of the system, lead- 
ing to an increase or a decrease in the amplitude of the 
oscillations. 



3. SAM or SLM? 

To get a strong SAM or SLM, a good beating contrast 
is needed. Then if the envelope has a minimum or a max- 
imum at the surface, it is a SAM or a SLM, respectively. 
Indeed, from Eq.fJTS]) the phase of the envelope relative 
to the SL surface depends essentially on <j) and on the 
sign of Q. The latter is seen in Fig[6]to be negative below 
the gap and positive above. In this section we focus on 
</>, the dephasing of the envelope function. 

For a given type of SL (materials and thicknesses), 4> 
depends essentially on the frequency v = ^ and on the 
reduced dephasing a of the SL. As in the whole article, we 
use here the numerical parameters of the SL described in 
Sect. Ill Dl whose dispersion diagram is plotted in Fig. [3| 

Figure [7K) presents in a 3D plot the variations of cj> as 
a function of a and of the frequency v, calculated using 
Eq. (JTTJ) around the second band gap. Figure[7p) reports 
the value of <f> as a function of a for two frequencies above 
and below the gap: 0.5807 THz and 0.6182 THz. These 
two frequencies reported in Fig. [6j correspond to a 0.5 
beating contrast C (Eq. (|B4])). so that at 0.5807 THz 



7 




0.5 3.55 0.5 0,65 0,7 

<J> V 




a 



FIG. 7: (color online) a) Value of <j> (color blue (f> — —n/2, 
red cj> = t/2) as a function of the frequency v (THz) and 
the reduced dephasing a. In the gap, <fr is not defined, b) 
(j> as a function of a at 0.5807 THz (black) and 0.6182 THz 
(red). Dark curves separate SAM and SLM regions. These 
curves have been calculated looking for couples (a, v) that 
equal the envelope wave function amplitude in z = and 
z = ^7, and correspond to <j> = ± 5 (cf Eq. fl!9|l ). Parameters 
corresponding to the SL described in section III Dl have been 
used. 



(resp. 0.6182 THz), the displacement u(z) (Eq. ([IB}) is 
dominated by the cos(Gz + ip/2) cos(kz + cf)) term (resp. 
sin(Gz + ip/2) sm(nz + cf)). 

From Fig.JTK) and b), 4> is a periodic function of a with 
period | at fixed frequency v. indeed, we are looking at 
the second band gap due to the harmonic at wavevector 
2G; by changing a from to 1, the surface z = sweeps 
two periods of that harmonic. For example, the two 
remarkable modes discussed above for which p2 is real, 
could be obtained for two values of a: these values arc 
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obtained by cancelling the d-i term in Eq. (|13c|) (n = 2). 
This has also been illustrated in FigO The case p 2 real 
and positive corresponds to a = a 2 = 3/4— 7/2 ~ 0.39 or 
a = «2 + l/2 = 5/4 — 7/2 ~ 0.89, whereas the case^2 real 
and negative corresponds to a = ai = 1/2 — 7/2 ~ 0.14 
or a = ai + 1/2 = 1 - 7/2 ~ 0.64. 

In Fig. [7^), in addition to the value of </>, we plot dark 
curves separating SAM and SLM regions. These curves 
have been calculated looking for couples (a, v) that equal 
the envelope wave function amplitude in z = and z — 
and correspond to (j) = ±j (cf Eq. (HU)). 

Below the gap, Figure [7£i) shows that regions of SAM 
dominate: especially, as the frequency closes with the 
gap, the SLM regions are reduced. This result is also 
coherent with Fig.[7j3): for 0.5807 THz, the slope of <fi as 
a function of a in the vicinity of a = 3/4 — 7/2 ~ 0.39 is 
very steep, and would be steeper for a frequency closer 
to the gap. Above the gap, the same conclusions apply: 
SAM regions dominate the phase space. 



4. SAM /SLM character 

In order to summarize the results of the two previous 
sections, let us define a "SAM/SLM character" as: 



X = In 



3 (0)/ Urms (|-) 



(20) 



which is positive for a SLM and negative for a SAM. The 
highest \x\, the more pronounced is the SAM or SLM 
character of the mode. \ is plotted as a function of the 
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FIG. 8: (color online) SAM/SLM character x as defined in 
Eq. (|20|l . A black curve delimits regions of positive and neg- 
ative x- Parameters corresponding to the SL described in 
section Hi Dl have been used. 

frequency v and the reduced dephasing a of the SL in 
FigEl It is clear that the more pronounced SAM or SLM 
(highest Ixl) are localized close to the gap, and that SAM 
dominate the phase space. 



5. Numerical Study 

To check our analytical predictions based on the SHA, 
we numerically integrate the differential equation Eq. ([3|) 
using a 7th-order Runge-Kutta integrator around the 2nd 
band gap: we thus take into account all harmonics of the 
inverse square sound speed. We choose here to represent 
two different cases closed to the remarkable modes we 
point out above: P2 real positive for a = «2 = 3/4 — 
7/2 w 0.39 andp2 real negative for a = a\ = 1/2 — 7/2 w 
0.14. In both cases, we compute solutions of Eq. for 
0.5807 THz and 0.6182 THz (0.5 beating contrast) with 
free boundary conditions at the surface, i.e. u(0) = 1 and 
35(0) = 0; these numerical results are plotted in Fig. [9] 
For p 2 real positive (a — a 2 ), we clearly observe a SAM 
above the gap and a SLM below, whereas for p 2 real 
negative (a — oti), it is the contrary. With respect to 
the SAM or SLM character, the numerical results thus 
fully agree with our analytical analysis. 

Fig. [9] not only validates our calculations and approx- 
imations, it also confirms that the parametric resonance 
analogy is fruitful to get the main physics of wave prop- 
agation near band gaps in semi-infinite SL. However, 
please note that numerical values of the contrasts C and 
the envelope wavelengths measured in Fig. [5] differ from 
the ones predicted by Eq. (|B4)) and Eq. (JTUJ. We at- 
tribute these discrepancies to the SHA. Indeed, the SHA 
allows the prediction of the width of the 2nd band gap 
with an error of the order of |pi| 2 . Thus, the gap edges 
predicted by the parametric oscillator analogy slightly 
differ from the one calculated from the exact transfer ma- 
trix method (or by a direct integration of Eq. ©). Due 
to the almost horizontal dispersion diagram curve near 
the gap edges, a small variation of frequency to may in- 
duce a non negligible variation in k. Consequently, while 
the parametric oscillator analogy gets the main physics 
and especially allows us to find out the parameters to see 
SLM, it remains based on the SHA and thus, some care 
should be taken in order to predict quantitative values, 
which will better be derived from an exact approach like 
the transfer matrix method. 



IV. OBSERVATION OF SLM AND SAM 

We now address the problem of the experimental ob- 
servability of NBPM. Trigo et alJ£ claim they have 
shown the existence of SAM using pump-probe experi- 
ments coupled to a theoretical analysis. In the following, 
we would like to focus in detail on the optical activity of 
NBPM. 

Colvard et a\J£ addressed the phonon-electron cou- 
pling in GaAs/AlAs SL. The reduced dephasing of 
phonons compared to layers of the SL is a key param- 
eter for the strength of this coupling. They showed in 
agreement with experiments that, out of electronic reso- 
nance conditions, B2 phonons, in phase with layers of the 
SL, do not create a strong electron-phonon coupling and 
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FIG. 9: Numerical solutions of Eq. © for 0.5807 THz (below 
the gap) and 0.6182 THz (above the gap) with free boundary 
conditions u(0) = 1 and ^f(0) = 0. The reduced dephasing 
a of the SL compared to the surface z — is respectively 
a = Q2 = 3/4 — 7/2 ~ 0.39 for figure a) and b) and a — 
ai — 1/2 — 7/2 ~ 0.14 for figure c) and d). For each graph, 
we plot both a zoom and overall picture of the mode. The 
overall picture exhibits the Loving or Avoiding character of 
the modes. The zoom allows the discussion on the phases 
of the mode compared to the SL when the amplitude of the 
mode is maximum. Dashed parts represent AlAs layers and 
blank ones GaAs layers. 



are thus not observable. On the contrary, Ai phonons, in 
quadrature phase with layers of the SL do create a strong 
electron-phonon coupling and are thus observable. These 
conclusions rely on the electron-phonon coupling mech- 
anism: the deformation potential (in GaAs/ AlAs SL), 
whose hamiltonian is proportional to the divergence of 
vibration modes. 

(2) 

In the case of NBPM around the band gap at w BG , 
their wave vector G±k is slightly different from the prim- 
itive vector G of the reciprocal lattice, thus their relative 
phase compared to layers of the SL slowly shifts with z 
as ±kz. Thus, such modes could be in phase or out of 
phase with the SL depending on the position z in the SL 
(as long as k ^ 0). However, though the envelope ampli- 
tude is varying, we can assume that these modes will be 
observable if they have the A\ symmetry (i.e. in quadra- 
ture phase with the layers of the SL) at their maximum 
envelope amplitude. 

As demonstrated in Appendix [C] for the remarkable 
cases of Sect. MI B 11 we find that, out of resonance con- 
ditions, observable modes are found below the gap inde- 
pendently on their SAM or SLM nature. 

To check our latter analytical prediction, numerical 
solutions of Eq. © below (0.5807 T Hz) an d above 
(0.6182 THz) the gap obtained in Sect. MI B 5l are com- 
pared to the SL in regions of maximum envelope am- 
plitude. Fig. [§] reports these results. In both cases 
P2 real positive (a = 3/4 — 7/2) and pi real negative 
(a = I/2—7/2), phonons above the gap are in phase with 
the layers of the SL at their maximum amplitude enve- 
lope in agreement with our prediction: they are hardly 
observable. On the contrary, phonons below the gap are 
in quadrature with the layers of the SL and will thus 
induce a strong electron-phonon coupling in agreement 
with our predictions. 



V. DISCUSSION 

We have shown that the envelope amplitude of near 
Brillouin zone center acoustical phonon in the vicinity of 
the surface can be maximum (SLM) or minimum (SAM) 
depending on the reduced dephasing of the SL relative 
to the surface. For a given function -^?(z), Eq. (fT|) 
is a second order differential equation whose solutions 
are entirely determined by two boundary conditions we 
can reduce to the definition of (ti(Q),4j(0)). Knowing 
one solution u±(z), U2(z) = Ui(z + z re f) is also a solu- 
tion for the boundary conditions 1*2(0) = u\(z re f) and 
*&(0) = ^(zref), and for the function £(2 + z ref ). 
Consequently, if SAM exists, a judicious choice of z re f 
allows to obtain SLM: one just needs to cut the SL at 
z re f chosen in the region where the envelope of the SAM 
is maximum. Assuming free boundary conditions, this 
choice imposes the reduced dephasing a of the SL, which 
precisely corresponds to our study. Reciprocally, choos- 
ing a = 0, the choice imposes the boundary conditions. 

We now address two discrepancies between our results 
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and related works in the literature. 

The first one concerns conclusions of Ref. [26]: "In the 
physically important cases of a free and a clamped in- 
terface, we emphasize, surface repulsion is a property of 
all modes with wavevectors close to those of dispersion 
gaps", suggesting that only SAM exist. We have shown 
in this work, that in the case of free boundary conditions, 
both SAM and SLM exist. The discrepancy between our 
results and the cited work may be explained since cal- 
culations of Ref. Hil only consider a SL with a reduced 
dephasing a — relative to the SL surface. In the pre- 
cise case of free boundary conditions and a = 0, our 
results agree with conclusions of Ref. |26|. 
The second discrepancy concerns claims of Ref. HH: 
NBPM "have a tendency to avoid the boundaries, ir- 
respective of the boundary conditions" suggesting that, 
whatever the boundary conditions, only SAM exist (note 
that this worki£ also only considers SL with a = 0). If 
we have restricted our study to the case of free boundary 
conditions, of course, the study could easily be extended 
to the case of any boundary conditions. Actually, as sug- 
gested in the begining of the section, we found that fixing 
the reduced dephasing a — and varying the boundary 
conditions produce the same qualitative conclusions as 
fixing the boundary conditions and varying the relative 
phase: both SAM and SLM exists and SAM regions dom- 
inate the phase space. We have chosen not to discuss 
in detail such effects of changing boundary conditions 
because the case of free boundary conditions is experi- 
mentally the most relevant. To illustrate the correctness 
of our analysis, and in disagreement with statements of 
Ref. [H, the reader can easily check the existence of a 
SLM by direct integration of Eq. © using the following 
parameters: u(0) = 1.5 and ^(0 ) = 1 and v = 0.5807 
THz and the SL defined in Sect. HlDl with a = (the 
same as Ref. 0) . 

Let's now consider experimental observations of SAM 
or SLM. Ref. [Til claims to observe SAM under the gap us- 
ing a SL with a = and free boundary conditions. This 
observation agrees with our analysis: from Fig. [5] NBPM 
under the gap for a = are SAM and from Sect. IIVI they 
are likely to be observable. 

In Sect. [TV] we show that both SAM and SLM could 
in principle be observable. We would like to address here 
technical problems associated with the observation of 
SLM. Indeed, in Fig. [SJ we show that SAM regions domi- 
nate the phase space. Thus, a high precision on the value 
of a is needed if one wants to observe SLM. This might 
prove challenging for SL growth. As an illustration, let's 
consider the precision on the thickness of layers produced 
by molecular beam epitaxy, which is about 2 monolay- 
ers («0.5 nm). In the case of the SL already studied 
(AlAs(2.35 nm)/GaAs(5.9 nm)) a precision of ±0.25 nm 
on the value of the dephasing r corresponds to a preci- 
sion of ±0.03 on the value of a. If a exactly matches 
3/4 - 7/2 = 0.39, the NBPM is a SLM for any value of v 
below the gap. However, if a lays in the range 0.39±0.03, 
The frequency v has to be lower than 0.583 THz to get 



SLM. Thus, because of the lack of precision on the layers 
thicknesses, observations of SLM would require a pro- 
tocol that selects NBPM with a low enough frequency 
(below the boundary SAM/SLM shown in Fig. [5]): we 
think such a selection may be achieved by the presence 
of a substrate. Indeed, the length of the SL and the pres- 
ence of a substrate is an important parameter that tunes 
the lifetime of observable phonons created in the SL by 
the excitation laser. A precise study of these effects will 
be detailed in a forthcoming publication. 

VI. CONCLUSION 

We have shown using a fruitful analogy with the para- 
metric oscillator that NBPM could either avoid (SAM) 
or love the surface (SLM) depending on the relative 
phase of the SL compared to the sample surface. More- 
over, we have shown that both modes should in principle 
be observable using judicious experimental parameters. 
Whereas SAM have already been observed, experimental 
evidences of SLM will be difficult in view of the techni- 
cal challenge to achieve a convenient SL. Despite these 
difficulties, we think such experiments are possible choos- 
ing an appropriate SL length. The experimental distinc- 
tion between SAM and SLM may be demonstrated us- 
ing pump-probe experiments with different penetrating 
length lasers. 

Finally, we would like to underline that our study of 
sound waves in SL could be generalised to any kind of 
waves. So that electroni c 27 ' 28 or electromagnetic*^ wave 
functions in electronic SL or Bragg mirrors also show a 
Surface Avoiding or Surface Loving character depending 
on the surface termination of the SL. 



APPENDIX A: RESOLUTION OF EQ.(6) IN A 
SEMI-INFINITE SL 

Finding the solutions (|9"aj) and (|9b[) amounts to diago- 
nalize the square matrix in Eq. ([5]). The eigenvalues are 
±in, defined by (fT0|) . Since we are focusing to the out- 
side of the gap, k is real and we may choose it positive 
(Changing n by — n is equivalent to permute constant 
A(B)+ and A(B)_ in Eq. ©): 

K= 2CV /J2 -N 2 ! (Al) 

recalling S — k ° from Eq. (fTTI) . \p2\ denotes the mag- 

nitude of p2- [pa| = y/p2p2- After diagonalization, the 
following relations are found between the constants in 
Eq. ©: 

B + = -C,e-'^A + (A2a) 
A- = (A2b) 

(A2c) 
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where ip is the argument of p2, so that P2 = \p2\ e 1 ^ and 



c 



s 

W 




(A3) 



Let us now impose a real amplitude Uq at the surface, 
u(0) = uq. This boundary condition reads 



A + (1 - (e^') + fl_ (1 - Ce +l,/ ') = «o 



(A4) 



We may try B_ = A + , since their respective factors are 
complex conjugates, and uq is real. Please note that 
changing uq to be complex would merely amount to de- 
phase the vibration in time. Then, introducing Aq (real) 

A e«2 +4 " 



and 



with 



> so that A + 

u{z) — Aq |cOS 

£cos 

^0 



, the displacement reads 



(G-K)z+$-(fi]\ 



Uq 



(cos (1 + 0) -Ccos (f-<£)) 



(A5) 



(A6) 



The boundary condition 4^(0) = yields to the fol- 
lowing expression which allows the determination of the 
integration constant 4>: 

/ , /„s (G + k)- C(G - k) 

tan cp = - tan (0/2 „, , ). „ ( A7 

(G + k) + C(G - k) 

Note that Eq. (IA7|) defines <f> modulo tt. However, u(z) is 
fully defined in Eq. (JTSJ) by defining the couple Aq and <j>. 
We thus use the convention to choose Aq > in Eq. (|A6|) : 
<j> is now defined modulo 2ir. Note that close to the gap, 
this procedure always yields a value of <f> in ] — % , 5]. 

We found one solution (|A5|) of this second order differ- 
ential equation. This is the only one satisfying the given 
free boundary conditions. Finally, Eq. (|A5|) can be writ- 
ten in the alternative form (|15|) , more convenient for the 
interpretation. 



APPENDIX B: ENVELOPE EQUATION 

As seen in FiglHl the displacement consist in a fast 
oscillating term at wave vector G, modulated by a lower 
wave vector envelope. One way to define this envelope is 
to take the square of u in Eq. (fT5l) 



l {z) = AAlUl - C) 2 cos 2 (Gz + |) cos 2 {kz + </>)+ 
(1 + C) 2 sin 2 (Gz + $) sm 2 (Kz + cf>)~ 
(1 - C)(l + C)2 cos(Gz + |) sin(Gz + f ) 
cos(kz + (f>) sin(«;z + 4>) j 



and then to average on the characteristic time of the fast 
oscillating term. So, we define: 



u rms( z ) = 2^ I - U 2 (z')dz' 



(B2) 



This procedure amounts to take a local mean, around 
every position z, on a ^ range. In Eq. (|B1|) , using k <C 
G, the cos 2 (Gz + f ) and sin 2 (Gz + f ) yield | and the 
last term vanishes. Finally, we end up with a local "root 
mean square" positive amplitude u rms 

u 2 ms (z) ~ 2A 2 [l + C 2 ~ 2C cos 2(kz + 0)] (B3) 

which will also be called the envelope of the displacement 
in this paper. The contrast of the beatings in u may be 
defined as 



C 



max(u rms ) - min(M rms ) 
max(u rms ) + min(w rms ) 



(B4) 



which is equal to £ above the gap and to ^ below the 

gap, since max(ii rms ) = ^2^0(1 + |C|) an( i niin(u rm s) = 
V%A \1~ |C||, with ( e]0, 1] above the gap and ( e] - 
00, —1] below the gap. 



APPENDIX C: OBSERVATION OF SAM OR SLM: 
ANALYTICAL STUDY 

To check if NBPM around band gap at u BG have the 
A\ symmetry at their maximum envelope amplitude, 
we examine their relative phase compared to the layers 
of the SL. So that, whereas in Sect. IIIIBI we were 
interested in the phase <j) of the envelope compared to 
the SL surface, we now get interested in the phase j of 
the quickly oscillating function (at spatial frequency G) 
compared to the layers of the SL. 



(Bl) 



Eq. (|T5|) gives the analytical expression of NBPM. The 
layers of the SL can be described using the inverse square 
sound speed as shown in Fig. [¥] However, only the first 
harmonic (at G) of the SL is necessary to get a pic- 
ture of the layers of the SL: the coefficient pi defined 
by Eq. (|13d|) determines the phase of layers of the SL. 

We thus, only have to compare the relative phase of 
the phonons (at their maximum envelope amplitude) to 
the one of the first harmonic of the inverse square sound 
speed. 

We now examine the two remarkable cases mentioned 
in Sect lBIBll 

Let us start with the cases P2 real and positive obtained 
for a = 3/4 - 7/2 or a = 5/4 - 7/2: the SAM mode, 
above the gap, varies like sin(Gz) sui(kz), whereas the 
SLM, below the gap, varies like cos(Gz) cos(kz). From 
Eq. (|13bp and Eq. (|13cp , we can see that c\ = and 
((i / which implies that the first harmonic of the SL 
varies like sin(Gz). Thus, it turns out that the SAM 
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mode is in phase with the SL at its maximum ampli- 
tude (kz = tt/2[tt])): the electron-phonon coupling is 
thus weak which makes it hardly observable. On the con- 
trary the SLM mode is in quadrature with the SL at its 
maximum amplitude (kz = 0[7r])): the electron-phonon 
coupling is high and so, this SLM mode is likely to be 
observable. 



A similar analysis for the cases P2 real and negative 
(a = 1/2 — 7/2 or a = 3/2 — 7/2) leads to the conclusion 
that in that case, the SLM mode, above the gap, will be 
hardly observable whereas the SAM mode just below the 
gap, is likely to be observable. 

Hence, in both remarkable cases, observable modes are 
found below the gap. 
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